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ABSTRACT 

Cosmic-ray energy densities in central regions of starburst galaxies, as inferred from 
radio and 7-ray measurements of, respectively, non-thermal synchrotron and 7r°-decay 
emission, are typically U p = 0(100) eVcm -3 , i.e. typically at least an order of magni- 
tude larger than near the Galactic center and in other non-very-actively star-forming 
galaxies. We hrst show that these very different energy- density levels reflect a similar 
disparity in the respective supernova (SN) rates in the two environments, which is not 
unexpected given the SN origin of (Galactic) energetic particles. As a consequence of 
this correspondence, we then demonstrate that there is partial quantitative evidence 
that the stellar initial mass function (IMF) in starburst nuclei has a low-mass trun- 
cation at ~2M , as predicted by theoretical models of turbulent media, in contrast 
with the much smaller value of 0.1M Q that characterizes the low- mass cutoff of the 
stellar IMF in 'normal' galactic environments. 
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In the nuclear regions of starburst galaxies active star for- 
mation powers emission of radiation directly by supernova 
(SN) explosions and indirectly by SN-shock heating of in- 
terstellar gas and dust, as well as from radiative processes 
involving SN-shock-accelerated electrons and protons. 

Timescales of starburst activity in galaxies are compa- 
rable to galactic dynamical timescales, tsb ~ Td yn ~ 10 8 yr. 
On the other hand, the timescales required for cosmic ray 
(CR) protons to gain energy by the Fermi-I acceleration 
process and to lose it (via pion decay into photons and 
e + e~ pairs, or advection) are r+ ~ r_ ~ 10 5 yr, respec- 
tively. Indeed, assuming the Fermi-I process to be at work 
in a SN remnant, r+ = E/E = (AE/E)~ 1 At = (3~ 1 At = 
(10//3o.i)Ai ~ 10 5 yr, where At ~ 10 4 yr is a typical SN 
remnant lifetime, and /3o.i is the speed of the SN ejecta in 
units of 0.1c. The main energy-loss timescale is largely de- 
termined by the starburst-driven outflow from the starburst 
region, i.e., typically r_ ~ T out ~ 10 5 yr (see Sect. 4). Thus, 

T + ~ T_ << T S B- 

The relative lengths of these timescales suggest that in a 
starburst region a balance can roughly be achieved between 
energy gains and losses for galactic CRs during a typical 
burst of star formation. Under basic hydrostatic and virial 
equilibrium conditions in a galaxy, a minimum-energy con- 
figuration of the field and the CRs may be attained [but see 
Beck & Krause (2005) for a critical view]. This implies that 



energy densities of particles and magnetic fields are roughly 
in equipartition (e.g., Longair 1994). 

The equipartition assumption enables deduction of the 
CR proton energy density, U p (essentially all the particle en- 
ergy content), from the measured synchrotron radio emission 
(which can be observed relatively easily) , and a theoretically 
expected ratio of the CR proton to electron energy density, 
thereby in principle circumventing the need for a second 
observable, such as the predicted - and only very recently 
detected (from the nearby starburst galaxies NGC 253 and 
M82: Abdo et al. 2009, Acciari et al. 2009, Acero et al. 2009) 
- very-high-energy (VHE: >100GeV) 7-ray emission. This 
emission is largely from CR proton interactions with am- 
bient gas protons, which produce neutral (ir°) and charged 
pions; 7r° decays into 7-rays. On the other hand, because in 
a starburst region r_ << tsb, U p can be roughly estimated 
also from SN rates and the fraction of SN energy that is 
channeled into particle acceleration. 

In this paper we argue the the basic underlying connec- 
tion between radiative yields of CR electrons and protons 
and SN rates, can be exploited to obtain some insight on 
the nature of the stellar initial mass function in the star- 
burst region. In Assuming particle- field equipartition we de- 
rive (Sect. 2) an approximate expression for U p as a function 
of the radio flux. Using radio data from the literature, in 
Sect. 3 we show that U p ~ 0(100) eV cm -3 in the central 
starburst nuclei of three well-known actively star forming 
galaxies. For the same starburst galaxies, in Sect. 4 we esti- 
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mate U p from the statistics of SN events and the (putative) 
SN energy that goes into particle acceleration. We show that 
these respective estimates are in reasonable agreement, even 
if appreciably uncertian. In Sect. 5 we outline how this con- 
sistency can be used as a diagnostic of the initial stellar 
mass function (IMF), particularly of its low mass end. Our 
conclusions are summarized in Sect. 6. 



2 PARTICLE AND MAGNETIC FIELD 
ENERGY DENSITIES 

Energetic electron populations consist of primary (directly 
accelerated) and secondary (produced via charged pion de- 
cays) electrons. Their combined spectral density distribution 
is assumed to be a single power-law of the Lorentz factor, 7, 
in some interval 71 < 7 < 72 



Table 1. Star-forming galaxies: IRAS data. 



iV e (7) = N e , (1 + x) r 



(1) 



where -/V e ,o is a normalization factor of the primary electrons, 
and x the secondary-to-primary electron ratio. 

Electron synchrotron emission from a region with radius 
r s (of an equivalent sphere that has the same volume) and 
mean magnetic field B, located at a distance d, is 

/„ = 5.67 x 1(T 22 ^JV e ,o(l + x) a(q) B^ x 



4 x 10 6 



q-i 



erg/(s cm Hz) 



(2) 



where a(q) is defined and tabulated in, e.g., Tucker (1975). 
Setting tp = (r s /0.1kpc) -3 (d/Mpc) 2 (/iGHz/Jy) and v = 
1 GHz, the the normalization of the electron spectrum 
(Eq.(2) is 

N efi = ^-xW- 15 tPa( q y 1 B- 3 i 1 250 1 ^- (3) 

A second relation is required to separately estimate N e 
and B. This is provided by equipartition between particles 
(electrons and protons) and the magnetic field, U p + U c ~ 
B 2 /8tt, which could possibly be attained due to a high de- 
gree of coupling between all the relevant degrees of freedom 
in the starburst region. In terms of the proton to electron 
energy density ratio, n(q)/(l + x) = U p /U c , equipartition 
translates to 



u P [i + 



1 + X ] 
K(q) J 



877 



(4) 



The electron energy density is U e = N e; o(l + x) rn c c 2 
■f-ri ^ 1 ~ <1 ^ n a PP rox i ma te expression for U c is obtained 
by ignoring the contribution of low energy electrons with 
7 < 71, whose energy losses are dominated by Coulomb in- 
teractions with gas particles. While the change of the elec- 
tron spectrum at low energies (when compared with that at 
higher energies, which is directly deduced from radio mea- 
surements) is important (e.g., Rephaeli 1979, Sarazin 1999), 
for our purposes here we can simplify by considering only 
the spectrum above 71. In the starburst regions of interest 
here 71 ~ 10 3 , but our numerical estimates depend only 
weakly on its exact value. Under this assumption, the elec- 
tron energy density is U e = iV e ,o (1 + x) rn c c 2 7 1 ~ 9 d7, 
where 72 is an upper cutoff whose exact value is irrelevant 



Object 


d <») 


J12/J. 


Ab) 
J 25m 


Ab) 
J 60m 


Ab) 
J 100m 


Ac) 

TIR 


NGC253 


3.0 


41.04 


154.67 


967.81 


1288.15 


44.00 


M82 


3.6 


79.43 


332.63 


1480.42 


1373.69 


44.35 


Arp 220 


72.3 


0.61 


8.00 


104.0 


112.0 


45.68 



( a ) Distances in Mpc. ( 6 ) IRAS flux densities (in Jy), from Gcnzcl 
ct al. (1998). M TIR luminosities (in erg s" 1 ), in log form. 

in the applicable limit 72 >> 71. For q > 2 and 72 >> 71, 
U e ~ iV e ,„(l + x)m oC 2 7 2 -V(g-2), which upon substitution 
of the expression for AT e>0 (l + x) from Eq.(3) yields 



U a = 



2.96 



x 1CT 22 250 5 tp ■ 



7i 



q+2 



B~ 



(5) 



(1 + X) T (q~2)a(q) 
Using Eq.(4), we then get the equipartition magnetic field 

250 9/2 ip 



7.44 x 10" 



7i 2 - 9 ' 



_ 

1 + X J (q-2)a(q) 



(6) 



Finally, using Eqs.(6),(4) and the expression for n(q) (e.g., 
Persic et al. 2008), we readily obtain an explicit expression 
for Ut>. 



2.1 Secondary electrons 

The ratio x is clearly a function of energy; it depends on 
the injection proton-to-electron number ratio, p/e, which 
generally depends on the injection slope (e.g., Schlickeiser 
2002), and on the gas optical thickness to pp interactions 
(that produce neutral and charged pions). Because of the 
high gas density in the starburst region the production of 
secondary electrons is expected to be relatively important. 

Each e + resulting from tt + decay (following a pp in- 
teraction) generates (through a photon pair obtained by in- 
teraction with an ambient e~) two e + e~ pairs. The e~'s 
from the pairs are secondary electrons whose energy is | 
that of the initial e + . The positrons in the pairs repeat the 
loop. This cascade continues as long as the photon energy 
is not smaller than 1.022 Me V. However, for magnetic fields 
~100/iG (see below) and an electron density of n c , this low 
threshold value is replaced by a higher threshold, E c , below 
which electrons lose energy much more rapidly by Coulomb 
interactions with gas particles than by synchrotron emission 
(e.g., Rephaeli 1979). With typical gas densities ~100cm~ 3 
in the central region of a starburst galaxy (such as M 82 and 
NGC 253, which we specifically consider in the next section), 
E c ~l GeV. In our approximate treatment here we use E c as 
the threshold energy for the production of secondary elec- 
trons. 

In order to evaluate the number of secondary electrons 
produced by a power-law distribution of CR protons where 
every proton undergoes one pp interaction, let us consider 
that a proton of kinetic energy T produces secondary elec- 
trons with energies successively lower by a factor 4, starting 
at T/3 and down to E c . 

In our simple shower model the number of secondary 
electrons produced per interacting proton of kinetic energy 
T - i.e. those produced down to the last step of the shower, 
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k max =[log(3S c T- 1 )/log(l/4)] + 1 - is 2 fcm "- 1 . (The first 
step of the shower, k=l, corresponds to the first e + result- 
ing from 7r + decay following the initial pp interaction.) In- 
tegrating over a oc T~ q CR proton spectrum and multi- 
plying by the relevant injection p/e ratio, we obtain the 
secondary-to-primary electron ratio corresponding to the 
case when all CR protons produce secondary electrons, 
Xo=r p/e fz in 2 [fcmax(T) - 11 T~' J dT, where r p/e is the proton 
to electron density ratio. For a strong non-relativistic shock 
(q=2, implying p/e=43), E c =l GeV, and _B in =10 TeV we get 
Xo— 50. However, not every CR proton produces secondary 
electrons. The mean free path of CR protons in a medium 
of density n p due to pp interactions is A pp = (fjp P rtp) _1 ; for 
protons with kinetic energy T ~ few TeV the cross section 
is <Tpp ~ 50 mb. For an ambient gas density n p ~ 100 cm~ 3 , 
App ~ 65 kpc. The probability for a single CR proton to 
undergo a pp interaction in its 3D random walk through 
a region with a typical radial extent r s ~l/4kpc is then 
V3r s /\ pp ~ 6.7 x 10~ 3 . We then estimate that for a strong 
non-relativistic shock the secondary to primary electron ra- 
tio is x=Xo V3 (r s /A pp )~0.33. 

In the rest of this paper we shall use x= 0-5 as a fidu- 
cial value. This is in approximate agreement with results of 
detailed numerical starburst models for energies <; 10 MeV 
(e.g., Rephaeli et al. 2009; see also Paglione et al. 1996; Tor- 
res 2004; Domingo-Santamana & Torres 2005; De Cea et al. 
2009). 



3 CR ENERGY DENSITY IN STAR-FORMING 
ENVIRONMENTS 

In the previous section we saw that if the source size, dis- 
tance, and radio flux and spectral index are known, then 
the corresponding particle energy density can be evaluated. 
As an application, in this section we compute the particle 
energy densities in the starburst nuclei of three nearby star- 
forming galaxies with different current star formation rates 
(SFRs) : the two nearby galaxies, NGC 253 and M 82, and the 
merging starburst Arp220 (SFR= 5, 10, and 220M Q yr~ 1 , 
respectively: see Table 1 and Eqs.(10),(ll)). These estimates 
are then compared with the corresponding mean value in the 
Galaxy [SFR~ (1 — 2) Mq yr -1 ]. 

For NGC 253 (distance d = 3 Mpc), following Rephaeli 
et al. (2009), we assume the following parameters for the 
central starburst emission: r B = 200 pc, /ighz = 5.6 Jy, and 
a ~ 0.75. Assuming that this radio flux is synchrotron emis- 
sion from a non-thermal population of energetic electrons, 
the spectral index is q~2.5, hence <z=0.0852 and k~8. Set- 
ting x=0.5, we obtain J7 p «75eVcm~ 3 . 

For M82 we adopt (see Persic ct al. 2008) d = 3.6 Mpc, 
r s = 260 pc, /ighz = 10 Jy, q~0.71, and q~2.42 (the latter 
implies a ~ 0.09 and k ~ 10). Setting x=0.5, we obtain 
[/p^97eVcm- 3 . 

In Arp220 (distance d=72 Mpc), the starburst activity 
takes place in a central molecular-gas disk of radius 480 pc 
and thickness 90pc that accounts for SFR~120 Mq yr -1 , 
and embeds (at a galactocentric distance of 200 pc) two ex- 
treme starburst nuclei of radii 68 and 110 pc that account 
for SFR= 50 and 35 Mq yr -1 , respectively. This scheme [de- 
rived from Torres (2004) and Eq.(10) below] is consistent 
with the galaxy- wide TIR-inferred SFR~ 225 Mq yr -1 . For 



use in Eq. (6) , for the disk alone we assume an effective spher- 
ical radius r s ~ 250 pc, /i ghz — 0.3 Jy, and a ~ 0.65: for 
X = 0.5 [see Torres (2004), the bottom of his Fig.7, for 
N p /N c = 100, consistent with q = 2.3 (Schlickeiser 2002) - 
the latter value implies a ~ 0.09 and k ~ 12], this leads to 
Up w 520 eV cm -3 . The two extreme starburst nuclei embed- 
ded in the disk are environments that are probably similar to 
very young SN remnants, and hence far from equipartition 
(e.g., Torres 2004). 

By contrast, in the Milky Way - a much more sedate 
star-forming environment with SFR~ (1 — 2) Mq - the 
value of the proton energy density is (7 p ~leVcm _3 at the 
Earth's position (see Schlickeiser 2002), and (6 ±3) eVcm~ 3 
in the ~200 pc region of the Galactic center, as inferred from 
the measured VHE 7-ray emission (Aharonian et al. 2006). 
In the even more sedate Large Magellanic Cloud (LMC), 
£7 P ~(0.2 — 0.3) eVcm -3 as inferred from the measured HE 
emission (Knodlseder 2009). 



4 COSMIC RAYS AND SUPERNOVA RATES 
4.1 Estimates 

Combining the SN frequency with the residence timescale of 
energetic protons that give rise to TeV emission in the star- 
forming region, and assuming a bona-fide value of the energy 
that goes into accelerating CR particles per SN event, we can 
obtain a second, independent estimate of U p . Moroever, as 
we show below this estimate for U p agrees well with that de- 
duced for each of the four galaxies in Sect. 3. This match does 
hint a link between cosmic rays and SNe, in agreement with 
previous suggestions (Ginzburg & Syrovatskii 1964; Salvati 
& Sacco 2008). 

The energy-loss timescale for pp interactions is r pp = 
(tTppCrip) -1 . For protons with kinetic energy T ~ few x 10 
TeV <7pp ~ 50 mb (hardly varying in the energy range of 
interest). Hence, 

T PP ~ 2 x 10 7 np 1 yr. (7) 

In the r s = 0.3 kpc central starburst of M 82 the gas mass is 
M ~ 2 x 10 8 Mq (Persic et al. 2008, and references therein), 
so n p ~ 100 cm" 3 ; therefore, r pp ~ 2 x 10 5 yr for O(10) 
TeV proton CRs in the central starburst of M 82. A similar 
estimate for NGC 253 [with r s = 0.2 kpc (Forbes & DePoy 
1992), Mh 2 = 7.2 x 1O 7 M (Mauersberger et al. 1996)] 
yields very similar value, r pp ~ 10 5 yr. In the star-forming 
disk of Arp 220, r pp ~ 1.7 x 10 4 yr (based on the gas density 
quoted in Table 2 of Torres 2004). 

However, for the galaxies considered in this paper the 
^-production timescale is longer than the CR proton out- 
flow timescale, v out . The latter is dominated by advection 
of the energetic particles out of the disk mid-plane region in 
a fast starburst-driven wind. For M82, v out ~ 2500 km s -1 
(Chevalier & Clegg 1985; Seaquist & Odegard 1991; Strick- 
land & Heckman 2009). Assuming a homogeneous distribu- 
tion of SNe within the starburst nucleus of radius r s , the 
outflow timescale is then 

r out = 3 X 10 4 ( 7 r^—) ( ^—r 1 yr. (8) 

V 0.3kpc' v 2500kms- 1 '' y V ' 

During the timescale in which protons actually reside 
in the central starburst region, namely the shorter of the 
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latter two timescales, r out , a number ^sNT ou t of SN explode 
and deposit the kinetic energy of their ejecta, _E G j = 10 51 erg 
(Woosley & Weaver 1995), into the interstellar medium. Ar- 
guments based on the CR energy budget in the Galaxy and 
SN statistics suggest that a fraction r;~0.05 — 0.1 of this en- 
ergy is available for accelerating particles (e.g., Higdon et 
al. 1998; Tatischeff 2008). We then express the CR proton 
energy density as: 

77 = «5 i/SN Tout 71 Ec > ( T ' s ~\ ~ 3 

P 0.3yr- x 3xl0 4 yr 0.05 10 51 erg ^0.3kpc' ) 

eV cm -3 . (9) 

From Eqs.(9),(16) we see that a substantial agreement 
with our equipartition estimates in Sect. 3 is reached for our 
five galaxies: 

NGC253: Taking r s =0.20kpc (Forbes & DePoy 1992) and 
^SN=0.12yr _1 (Rieke et al. 1988)), we get U p ~75 eV cm -3 ; 

M82: Taking r s =0.26kpc (Persic et al. 2008 and refer- 
ences therein) and ^sN=0.25yr _1 (Rieke et al. 1980), we 
get L/p~95eVcm~ 3 ; 

Arp220: Taking r s =0.25kpc (see Sect. 3) in Eq.(9) and 
0.09 kpc (i.e., the thickness of the starburst disk, perpen- 
dicular to which the wind will break out) in Eq.(8) and 
j/ SN «3.5yr~ (Smith et al. 1998; Lonsdale et al. 2006), we 
get ?7p~505eVcm~ 3 . 

Milky Way: for the relatively quiet star forming environ- 
ment of the Galaxy (SFR~ (1 — 2) M Q yr" 1 ) we consider 
a nuclear star-formation volume comparable to those con- 
sidered for the previous galaxies, r B = 0.2 kpc. The large- 
scale distribution of the HI gas in the Galaxy is made up 
of two major components (Kalberla & Dedes 2008): most 
prominent is an exponential HI disk with Rd = 3.75 kpc, 
coplanar to the stellar disk and that can be traced out to 
35 kpc. Surrounding this region there is a patchy distribu- 
tion of highly turbulent gas that reaches higher scale-heights 
and larger radial distances. If bhi = 1 cm -3 in the solar 
neighborhood, i.e. at a Galactocentric distance Ro = 8 kpc 
(Dehnen & Binney 1998), the average value within r s is 
nui — 8 cm -3 . Hence the relevant timescale to be used 
in Eq.(9) is r pp = 2.5 x 10 6 yr. The full disk Galactic 
SN rate, z/ SN = 0.02 yr -1 (Diehl et al. 2006), should be 
rescaled to the r B region: assuming an exponential stellar 
disk with R d = 2.5 kpc, (Binney et al. 1997 from COBE 
L-band data; see also Dehnen & Binney 1998), the mass 
contained within r B is 3.66 x 10 of that contained within 
Ro, hence the resulting SN rate to be used in Eq.(9) is 
0.73 x lO^yr" 1 . Assuming rj = 0.05, from Eq.(9) we fi- 
nally obtain U p ~ 5.7eVcm~ 3 . This value is fully consistent 
with that derived, for the same region, by Aharonian et al. 
(2006) from the diffuse VHE 7-ray emission measured there 
with the H.E.S.S. Cherenkov array. 

LMC: This satellite of the Milky Way has a very quiet star 
forming environment, characterized by vsn = 0.002 yr -1 
(van den Bergh &: Tammann 1991) in the star forming re- 
gion of the galaxy. The latter can be modeled as a trun- 
cated disk/spheroid with r t ^3kpc whose thickness is also 
~3kpc (Weinberg & Nikolaev 2001), we use r s = 3 kpc 
in Eq.(9). There is no mass outflow from the LMC, hence 
the energy-loss timescale for pp interaction is relevant here; 
with an average gas density of n p ^2cm~ 3 , this timescale 



is T_=T pp ~10 7 yr. Using these values in Eq.(9), we obtain 
C/ p «0.2eVcm . This value is consistent with that inferred 
from Fermi/LAT measurements of the average emissivity 
spectrum of the LMC (either including or excluding the 
bright star- forming region of 30Doradus), which suggest a 
Up value about 0.2-0.3 times that in the solar neighbourhood 
(Knodlseder 2009). 

In summary, the CR energy densities estimated in five 
galactic nuclei of similar size (three starburst galaxies, the 
central Galactic region, and the LMC), appear to be largely 
correlated with key features of the ongoing star formation: 
local ambient density, SN rate, and CR proton residence 
time. In fact, for same J?-E e j (by assumption) and similar r s 
(from observations), the CR energy density U p seems to be 
well described just as a function of the number of SN explo- 
sions during the CR residence timescale, ^sr-jTrcs , where the 
CR residence (energy-loss) timescale t^ b is defined in terms 
of the geometric mean of the pp interaction timescale and of 
the advection timescale, t^ s = t~ p (nm) + r~ ut (r s ). 

4.2 Uncertainties 

Uncertainties in the quantities appearing in Eq.(9) need to 
be assessed before we make deduction about the stellar IMF. 

• The values of U p derived in Section 3, based on particle- 
field equipartition and representing mean values for the re- 
spective central starburst regions, are essentially confirmed 
by accurate numerical treatments based on the solution 
of the diffusion-loss equation for the accelerated particles 
(M82: Persic et al. 2008, and De Cea et al. 2009; NGC253: 
Paglione et al. 1996, Domingo-Santamarfa & Torres 2005, 
and Rephaeli et al. 2009; Arp220: Torres 2004). Starting 
with the initial accelerated spectrum in the central starburst 
region, these authors followed the evolution of the spectrum 
as the particles diffuse and are convected from this inner re- 
gion to the outer disk (and halo). Including all the relevant 
leptonic and hadronic interactions in a numerical treatment, 
the particle spatial and energy distributions, and their re- 
spective radiative yields, were followed across the galaxy 1 . 

At the time of this writing 7-ray detections of M 82 and 
NGC 253 in the HE band (Abdo et al. 2009) and in the VHE 
band (Acciari et al. 2009; Acero et al. 2009) were announced. 
The measured fluxes and spectra of both galaxies in the two 
bands agree with the predictions of recent numerical models 
(Persic et al. 2008; de Cea Del Pozo et al. 2009; Paglione 
et al. 1996; Domingo-Santamarfa & Torres 2005; Rephaeli 
et al. 2009), which are all based on values of U p that are in 
agreement, to within a factor of w2, with our estimates in 
Sect.4. 

Given these arguments, we feel confident that the en- 
ergy densities of CR protons in starburst nuclei, as de- 
rived from our equipartition arguments above, are reliable 
to within a factor of J; 2. 

1 Only in the two extreme starburst nuclei of Arp 220 - which we 
don't consider in this paper - is particles/field equipartition most 
likely not attained. These strongly magnetized regions (where 
typically the mean field is at the mG level) are similar to very 
young SN remnants, and hence perhaps still too 'hot' for any sort 
of equipartition to be attained (Torres 2004). 
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• SN rates in starburst nuclei are in principle observable 
quantities. However, heavy optical extinction drastically re- 
duces one's ability to track SN activity, whereas radio counts 
of SN remnants need information on their ages in order to 
be turned into actual SN rates. In general, long monitoring 
timescales, with measurements of actual source fadings and 
new appearances, are required to estimate the actual SN rate 
from the observation of the compact radio sources (Ulvestad 
& Antonucci 1997). As for our three starburst galaxies, we 
assess (based on published observational results) that their 
SN rates are known to within a factor of ^5 1.5. 

• The uncertainty in the CR proton outflow timescale, r ou t, 
arises from the uncertainty in the fast wind velocity, v out . 
The latter is probably known to within ~50%. Indications 
of such velocity in M 82 come from multifrquency radio mea- 
surements of its nonthermal halo emission (Seaquist & Ode- 
gard 1991). The halo emission is interpreted as synchrotron 
radiation from relativistic electrons accelerated in SN rem- 
nants and then swept out the disk plane by the extensive 
wind associated with the starburst (Chevalier & Clegg 1985; 
Seaquist & Odegard 1991). The steepening of the radio spec- 
tral index (from — 0.4±0.1 in the starburst nucleus to about 
— 1 in the halo), is explained by a simple model which in- 
volves outward convection of relativistic electrons which suf- 
fer energy losses by Compton scattering against IR photons 
(which are plentiful in the presence of the central starburst) 
and adiabatic expansion (which makes this loss term largely 
dominant in the outer radio halo). In the context of this 
model a convection speed v out ~ (1000 — 3000) kms -1 is 
deduced, as well as a fraction of the total energy input to 
the halo in the form of relativistic particles of a few percent, 
consistent with the fraction found in Galactic SN remnants - 
this agreement, too, supporting the hypothesis that the rel- 
ativistic particles in the halo are produced by SN remnants 
in the disk (Seaquist & Odegard 1991). The relationship be- 
tween this fast particle outflow and the slower, denser minor- 
axis outflow observed at optical and X-ray wavelengths (with 
deprojected velocities of ~600kms _1 ) is unclear. Much of 
the optical data may be explained from large-scale shocks 
from the high-speed wind plowing into the gaseous halo and 
entrained disk gas (Shopbell & Bland-Hawthorn 1998). 
The assumed wind speeds are in agreement with the sug- 
gestion, from standard wind theory, that wind terminal ve- 
locities usually amount to some (small) multiple of the es- 
cape speed, v CBC (e.g., Parker 1958). For example, the ter- 
minal wind velocity may be zz(3 — 5) v csc in M82 (Strick- 
land & Heckman 2009), and «3w OS c in our Galaxy (Everett 
et al. 2008). Would the correspondingly short convective 
timescales imply similarly short r s lengthscales? For the real 
galaxies considered in this paper, r s are deduced from high- 
resolution optical and radio data, so their values are quite 
accurate. For distant galaxies, where such measurements are 
more difficult, the error must be evaluated theoretically and, 
conservatively, it may be a factor of ,$ 2. 

• Published estimates of the energy of a core-collapse SN 
available for Galactic CR acceleration, (0.5—1) 10 50 erg, are 
fairly robust - within a factor i; 2 (Lingenfelter et al. 1998; 
Higdon et al. 1998; Tatischeff 2008). This is ~5% - 10% of 
the total kinetic energy of the SN ejecta, ~10 jl erg (Woosley 
& Weaver 1995), that is also a commonly assumed nominal 
value for this parameter in all core-collapse SNe. We assume 



this 'Milky Way normalization' to hold also for the starburst 
galaxies considered in this paper. The observational fact that 
the much more massive (> 25 M ) Wolf-Rayet stars seem 
to explode with _E e j~10 52 ergs (e.g., Lozinskaya & Moiseev 
2007) is a subtlety that, although important, does not sub- 
stantially alter the general picture given the relative rarity 
of these stars. 

We conclude that the quantities of Table 1 galaxies rel- 
evant to Eq.(9) are quoted in the literature as observa- 
tionally precise, to within a factor of i; 2. In this case, 
the discrepancy between the very high CR energy densities 
deduced for active starburst nuclei ([/p«100eVcm~ 3 ) and 
the low value deduced for the nuclear region of the Milky 
Way ([/p«5eVcm ) seems very significant. Indeed, within 
our sample this discrepancy is especially striking between 
NGC 253 and the Milky Way, which share very similar TIR 
luminosities and thus star formation rates. 



5 SUPERNOVA RATES IN STAR-FORMING 
ENVIRONMENTS 

In this section we investigate the constraints on the stellar 
initial mass function (IMF) that could possibly be derived 
from the SN rates that have been observed or deduced for 
our sample galaxies. 

A galactic SN rate can be modelled once we know i) 
the total gas mass that is turned into stars each year (i.e., 
the SFR), and it) the fraction of young stars that have M > 
8Mq and are therefore SN progenitors. 

To evaluate the SFR we use Kennicutt's (1998a) popu- 
lar I/tir. — > SFR conversion: 

-t/TIR 



SFR 



2.2 x 10 43 erg 



where Ltir the total IR [i.e., ( 
turn, the TIR flux is given by 



Meyr- 1 , (10) 
- 1000)/im] luminosity. In 



/tir = 1.8 x 10" 11 [13.48/ia + 5.16 / 2 5 + 2.58 / 60 + /ioo] 

(11) 



-2 -1 

erg cm s 



(e.g., Sanders & Mirabel 1996), where fx are the IRAS 
spectral densities (in Jy) at A (in /im). The continuous- 
burst model of (10-100) Myr duration of Leitherer & Heck- 
man (1995) and the Salpeter stellar IMF with mass lim- 
its (0.1-100) Mq, provides a Ltir — > SFR conversion that 
applies most exactly to young starbursts embedded in op- 
tically thick dust clouds, such as the starburst regions in 
our sample galaxies. The calibration of Kennicutt's con- 
version has a ~30% uncertainty that lies within the range 
of published calibrations (see Kennicutt 1998b) 2 . In Ta- 
ble 1 we report all relevant values for NGC 253, M 82, and 
Arp220; the corresponding SFRs are, respectively, 4.5, 10.2, 

2 Observationally, the instantaneous IMF in a starburst is pos- 
sibly different from the time-integrated IMF (see review by 
Elmegreen 2005). Basically, because of the proto-stars formed in 
an instantaneous burst of star formation, the low-mass ones take 
a very long time to get to the main sequence, and by the time 
they do, the upper-main-sequence stars have already left the main 
sequence. This problem, of varying importance according to the 
type of environment being considered, can affect the calibration 
of Eq.(10) (see discussion in Kennicutt 1998b). 
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and 220 Mq yr . By comparison, our Galaxy has Ltir ~ 
iiboi ~ 5 x lO^ergs" 1 and SFR^ 2M yr^ 1 . 

If the SN progenitors are all stars with M>8M , the 
gas mass that is turned into SN progenitors is 



M gas (> 8M ) = SFR x 



8M 



- 1 dm 



m -x+i dm 



(12) 



where Mi ow and M up are, respectively, the lower and upper 
mass cutoffs of the stellar mass function, assumed here to 
be power-law with index x. 

Mass conservation requires that the gas mass turned 
into SN progenitors per year must equal the mass contained 
in >8Mq stars formed per year, 



M*(> 8M ) = A 



f 



dm 



(13) 



where A is the normalization of the mass spectrum of the 
stars formed per year. By equating the two mass rates we 
can determine the value of A. We can then compute the 
number of SN progenitors that are formed per year, which 
in steady state is equal to the number of SN explosions per 
year, dN S n/dt = A J^ p m~ x dm, that is 



z^sn = 



dt 



8M 



m x dm 



2.2 x 10 43 



dm 



(14) 



In the evaluation of Eq.(14) we assume a Salpeter (1955) 
stellar IMF, i.e. x — 2.35. The standard stellar mass range 
is usually taken to be bracketed by JVfi ow ~ 0.1 M and 
M up ~ 100 M Q (e.g., Tinsley 1980). The two mass limits 
are dictated by considerations of ignition and sustainabil- 
ity of nuclear reactions (Mi ow ) and of stability to radiation 
pressure (M up ). In what follows the exact value of the latter 
is relatively unimportant, whereas the value of the former 
is of great importance and hence will be discussed in some 
detail with reference to the physical status of the medium 
in the local star- forming environment. 

In a non-turbulent environment, arguments of grav- 
itational instability, efficient heat dissipation on collapse 
timescale, and thermal radiation of the collapsing object 
lead to the fragmentation limit, 

1/4 / N 1/4 

/ he 



M frag = 6/i. m p 



y jj,rn v c 2 I y Grrip 



(15) 



(e.g., Peacock 1999). For T ~ 100 A", it is M frag = 0.01 M Q , 
almost an order of magnitude lower than the hydrogen- 
burning limit. (Brown dwarfs and planets may form with 

Mf rag < M < Mburn, of COUrSe.) 

In a starburst environment, the (molecular gas) tem- 
perature is higher and the turbulence stronger. Because of 
the T 1 / 4 behavior in Eq.(15), the higher temperature causes 
only a modest increase of Mf rag . However, the substantial 
turbulence (comparable in energy to the gravitational en- 
ergy) in the highly inhomogenous environment, where all the 
stars form, explode, drive shocks, accelerate cosmic rays and 
drive a dynamo (McKee & Ostriker 1977), causes dense cores 
to be formed by turbulent shocks, and core sizes to scale 
with the thickness of post-shock gas, rather than with the lo- 
cal Jeans length. Therefore, fragmenting cores, much larger 
than the local fragmentation limit, can be assembled dynam- 



ically; indeed, the mass spectrum of dense proto-stellar cores 
generated by turbulent fragmentation is Salpeter-like for 
M> (1 - 2) M Q (Padoan & Nordlund 2002). The case for a 
top-heavy stellar IMF in a warm, dense, turbulent medium, 
due to a larger Jean mass, was also made by Klessen et al. 
(2007). 

Therefore in evaluating Eq.(14) we assume M up = 
100 Mq and Mi ow = 2 Mq . These assumptions yield 



0.12 yr 
0.26 yr 
5.6 yr~ 



NGC 253 
M82 
Arp 220 



(16) 



(For Miow = 3M a , the estimated usn values are, respec- 
tively, 0.14, 0.32, 6.8yr _1 .) These estimates are in excellent 
agreement with observational values: vsn — (0.1 — 0.2) yr -1 
(Antonucci & Ulvestad 1988; Rieke et al. 1980, 1988), and 
certainly usk < 0.3 yr -1 (Ulvestad & Antonucci 1997), for 
NGC 253; z/ S n ^ (0.2 - 0.3) yr^ 1 (Kronberg et al. 1985; 
Rieke et al. 1980) for M82; and u S n ~ (4 ± 2)yr _1 (Lons- 
dale et al. 2006) for Arp 220. Using the 'standard' low value 
Miow = 0.1 M Q , we would instead get u S k= 0.03, 0.08, 1.6 
yr -1 , which seems to be inconsistent with observations. 

The above result suggests that in the central starbursts 
of M82, NGC 253 and Arp 220, the observed SN rates are 
consistent with the respective TIR luminosities (which are 
proxies for their SFRs), if their stellar IMFs (assumed to 
be Salpeter-like) are "top-heavy" with the lower mass cut- 
off of Miow ~ (2 — 3)Mq. This deduction rings in with 
earlier suggestions on a high Mi ow truncation (at ~3M ) 
of the IMF in M82's nuclear starburst, based on the col- 
ors and luminosity of the galaxy (Rieke et al. 1980, 1993; 
Doane & Mathews 1993). Top-heavy IMFs have been ob- 
servationally suggested for Arp 220 (Parra et al. 2007) and 
IC 694 (Perez- Torres et al. 2009) based on the unusually 
large fraction of core-collapse SNe observed in these galax- 
ies. In comparison, in the Galaxy the 'standard' low-mass 
cutoff Miow = 0.1 Mq applies. With this low mass cutoff and 
_Ltir ~ 5 x 10 43 crgs _1 in Eq.(10), the resulting predicted 
SN rate is ^sn ~ 0.015 yr -1 , in agreement with the observa- 
tional value of ^sn ^ (0.02 ± 0.01) yr" 1 (Diehl et al. 2006). 
In the LMC, too, the 'standard' low-mass cutoff seems to 
apply; using Mi ow = 0.1 M and SFR«0.25 M yr _1 (from 
TIR, thermal radio, and Ha luminosities; e.g., Hughes et al. 
2007 and Whitney et al. 2008), the predicted SN rate in the 
LMC is 2-'sN~0.002yr~ 1 , in agreement with the observation- 
ally deduced value (van den Bergh & Tammann 1991). 

An immediate consequence of a top-heavy stellar IMF 
is that the SFR as deduced from the TIR luminosity ac- 
cording to Eq.(10) is overestimated with respect to its 
real value. The degree of overestimation depends on the 
detailed shape of the actual stellar IMF. As an illus- 
tration, consider a stellar IMF denned as </>(m)=0 for 
0.1M Q <m<Mi ow and 0(m)cxm~ 2 - 35 (i.e., Salpeter-like) for 
Mi ow < m< 100 Mq, then the real SFR would be only a frac- 
tion t'^° M@ m~ x+1 dm/ L 1 *? ,^ m~ x+1 dm of the nominal 

JM low ' J 0.1 Mq 

(i.e., Kennicutt-based) SFR: e.g., Mi ow =2M would en- 
tail SFR(real) = 0.29 x SFR(nominal). On the other hand, 
if more realistically (e.g., Padoan & Nordlund 2002) the 
IMF has a Salpeter-like profile for Mi ow >m>100M Q but 
it turns over at m~Mi ow such that <f>(m)ctm OL (with 
q>0) for 0.1 M Q <m<Miow, then the real SFR would only 



Probing star formation with galactic cosmic rays 7 



be a fraction f=[f™ m'-dm/ f™°™° m'-'dm] + 

M^/Cm" - 1 -^] / X m ° +ldm ^ the n ° minal 
SFR; with Mi ow =2M it is /=0.35, 0.36, 0.37, 0.42 for 
a=0,0.5, 1,2. 

The realization that the SFR deduced from Eq.(10) may 
be overstimated does not, however, invalidate our estima- 
tion of the SN rates in starburst environments. Underlying 
Eq.(10) is the assumption of a Salpeter IMF for the stars 
in the mass range relevant to heating dust clouds (i.e., stars 
with m ~ 5 Mq which emit copious UV radiation, where the 
absorption cross-section of dust grains peaks), which pro- 
vides the correct normalization of the massive star formation 
rate. Only as a further step an extrapolation down to 0.1 Mq 
is made in order to obtain the nominal SFR(>O.1M0) given 
by Eq.(10). Using a different calibration, say SFR(>5M ) 
as is common practice with radio SFR indicators (e.g., Con- 
don 1992), this numerical issue would not show up. How- 
ever, a revision - by a factor of ~3 — 4 down - of the SFR 
in the actively star forming galaxies at redshifts 1, when 
the cosmic SFR was much higher than the current average, 
would have a major impact on our quantitative understand- 
ing of the star formation history of the Universe (e.g., van 
Dokkum 2008). 

In summary, assuming the stellar IMF to be univer- 
sal (Salpeter-like in the present calculation), wc can model 
the fraction of newly formed SN progenitors (M > 8M a ) 
by varying the low-mass cutoff from the 'canonical' value of 
0.1 Mq, appropriate for a quiet environment like our Galaxy, 
to the 'top-heavy' value of (2 — 3) Mq, suggested for a turbu- 
lent starburst environment. For a given gas mass turned into 
stars per year in a galaxy, the fraction of SN progenitors can 
then vary by a factor of ~4. This variation is quite substan- 
tial, apparently well above the uncertainty in the quantities 
that enter Eq.(9). 



6 CONCLUSION 

The idea of a link between star formation and CR particles 
has been around for a long time. It is based on the recogni- 
tion that the CR electrons, produced in the sites of SNe ex- 
plosions (marking the end of >SMq stars progenitor stars), 
were responsible for diffuse synchrotron emission. The rough 
agreement between the relatively short lifetime (J; 3x 10 7 yr) 
of these massive stars, and the similarly short synchrotron 
energy loss time of the electrons, led to the expectation that 
the non-thermal radio emission of a galaxy is a measure of its 
star formation activity on scales much shorter than the Hub- 
ble time. Biermann (1976) first included the non-thermal 
radio emission, Lr, in quantitative galaxy population mod- 
els. A Lr-i/sn relation was given, theoretically, by Ulvestad 
(1982) and later, empirically, by Condon & Yin (1990). A 
starburst model, that unified features of the earlier models 
and used Lfir as the best estimate of the bolometric lumi- 
nosity of a starburst, was proposed by Gehrz et al. (1983). 
Finally Condon (1992), relating Lr and Lfir by attribut- 
ing both luminosities to massive (>5Mq) stars, proposed 
a widely used, although semiempirical (e.g., no mechanism 
for CR generation in SN remnants, nor CR propagation, are 
explained), radio SFR indicator (see Persic & Rephaeli 2007 
and references therein for further developments). 



The tight FIR/radio correlation (e.g., Condon 1992 and 
references therein) has been interpreted as a manifestation 
of the (apparently universal) link between non-thermal radio 
emission and recent star formation in galaxies. In particu- 
lar, if the dust contributing to Lfir is primarily heated by 
>5Mq stars (Devereux & Young 1990), if the CR energy 
production per SN is the same in all galaxies - starburst 
and normal alike - (Volk et al. 1989), and if the magnetic 
field strength either plays a minor role or it varies little be- 
tween galaxies (whether ordinary disk galaxies or luminous 
nuclear starbursts), then both Lfir and the CR production 
rate are directly proportional to the recent massive SFR. 

Finally, the use of Lfir (which in this paper is denoted 
by Lir) as an indicator of the current SFR has been long de- 
bated. The conversion is derived using synthetic population 
models. In the optically thick limit, it is only necessary to 
model the bolometric luminosity of the stellar population. 
In general, the choice of a suitable timescale for the stellar 
population (age of the starburst, timescale for sputtering of 
dust) will affect the calibration of the conversion: different 
assumptions by different authors have led to slightly differ- 
ent calibrations (see discussion by Kennicutt 1998b). For 
example, the conversion in Eq.(10) applies only to young, 
optically thick starbursts; whereas in more quiescent nor- 
mal star-forming galaxies the relation is more complicated, 
due to the increasing contribution to dust heating from from 
old stars, and the decreasing dust optical depth. 

In the present paper we have explored the link between 
star formation and CR protons; specifically, the link be- 
tween ongoing star formation and Lp. Doing this we have 
followed up, and to some extent ultimately strengthened, 
the longtime suggestion that CRs are related to SN explo- 
sions (Ginzburg & Syrovatskii 1964; see also Salvati & Sacco 
2008). 

We started by assuming that equipartition arguments 
provide realistic estimates of the particle energy densi- 
ties, Up, in some well-known starbursts with measured syn- 
chrotron radio emission. We then found that in these envi- 
ronments, for a given gas mass per year turned into stars 
(measured by Ltir), a "top-heavy" stellar IMF with a low- 
mass cutoff of ~2 Mq (as befits protostar fragmentation 
in a turbulent medium) leads to the nominal SN rate re- 
quired to sustain the correct U p . By contrast in our Galaxy 
and in the LMC, where star formation is much more mod- 
erate, the (measured) U p is consistent with the standard 
(Miow = 0.1 M©) IMF. In general, we may conclude that 
the instantaneous IMF in a starburst environment is possi- 
bly different from the time-integrated IMF (the latter being 
represented by the standard Salpeter function). 

A top-heavy stellar IMF should have appreciable effects 
on galaxy colours. For a given M gas , hot massive OB stars 
are more numerous, so the IR hump in the galaxy's broad- 
band continuum is enhanced, if the IMF is top-heavy rather 
than standard. Starburst galaxies are then expected to turn 
up more IR-bright than they would be with a standard IMF. 
Indeed, a detailed analysis of the colors of M 82 and NGC 253 
led Rieke et al. (1980, 1993) and Doane & Mathews (1993) to 
suggest that the stellar IMF there was truncated at Mi ow ~ 
3M . 

In each starburst galaxy considered in this paper, for the 
relevant quantities appearing in Eq.(9), we have used nom- 
inal values from the literature. These, although statistically 
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consistent within a factor of J; 2, may however be affected 
by systematics that are presently difficult to model, so the 
quantities appearing in Eq.(9) could be known more poorly 
than we may think. Given the only initial direct 7-ray-based 
determinations of U p , and the errors (statistical and system- 
atic) affecting the relevant quantities (z/sn, t_, rjE c - s , r s ), our 
present deductions on the stellar IMF are clearly only ten- 
tative. 

Indeed the question whether a low mass limit on the 
stellar IMF in a starburst region can be derived has been 
argued to be ultimately difficult (e.g., Elmegreen 2005), be- 
cause we are probably missing some crucial subtlety in the 
time-integrated star formation. Perhaps a relatively robust 
result is that, due to the extra energy input available from 
the starburst environment, the timescale required for proto- 
stellar clouds to collapse, get rid of their angular momentum, 
split up, and merge again, is just very much longer within 
a starburst than outside the region. So the corresponding 
low-mass stars would presumably stay in the Hayashi phase 
for a time longer than the starburst, and would arrive on the 
main sequence only after the starburst activity has subsided. 

In any case, as concerns the derivation of the massive 
SFR from CR protons, direct 7-ray detections of starburst 
galaxies (in the 0.1-100 GeV photon energy range with the 
Large Area Telescope onboard the Fermi Gamma-ray Space 
Telescope; and in the >100 GeV range with Cherenkov tele- 
scope arrays such as H.E.S.S., VERITAS, and MAGIC) have 
just started providing direct U p measurements, which hap- 
pen to be consistent with theoretical predictions, for at least 
some of the nearby galaxies. Given these initial agreements, 
should our predictions on U p be more generally confirmed, 
some immediate implications would be: 

(i) star forming galaxies are powerful particle accelerators, 
able to achieve CR energy densities orders of magnitude 
higher than what is measured in the Milky Way and in the 
LMC; 

(ii) SNe, both in quietly star- forming galaxies like the Milky 
Way and the LMC and in very actively star-forming galaxies 
such as those considered in this paper, probably share a 
common (universal?) CR acceleration efficiency; 

(Hi) particle energy densities and equipartition magnetic 
fields can be used (at least statistically) as proxies of the 
actual quantities; this could be particularly useful in the 
case of galaxies that are too far for their (unbeamed) 7-ray 
emission to be measured; 

(iv) the initial stellar mass spectrum shows a systematic 
dependence on the local star-forming environment (see also, 
e.g., Elmegreen 2005); and 

(v) contrary to common practice, each given SFR indica- 
tor should have different calibrations according to whether 
galaxy luminosities are dominated by the current or the av- 
erage past SFR - the former case including local starburst 
galaxies, ultra-luminous IR galaxies (ULIRGs), and high- 
redshift galaxies; the latter case including most local galax- 
ies (e.g., Persic & Rephaeli 2007). If so, since in the earlier 
Universe star formation took place in denser environments 
than at the present epoch, the calibration of SFR indica- 
tors would show a systematic dependence on redshift. This 
would deeply affect our understanding of the star formation 
history of the Universe (van Dokkum 2008) . 
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